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Abstract 



Slavnov-Taylor is a Mathematica package which allows us to perform automatic 
symbolic computation in BRST formalism. This article serves as a self-contained 
OO \ guide to prospective users, and indicates the conventions and approximations used. 
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CLi' 1 Program summary 

Title of program: Slavnov-Taylor 
^ \ Available at: http://pcpicariello.mi.infn.it/ST/ 

' Programming Language: Mathematica 4.0 

Platform: Any platform supporting Mathematica 4.0 
Computers tested on: Pentium PC 

Operating systems under which the program has been tested: Linux 
Memory required to execute: Minimal: 1.254.784 bytes, Standard: 1.281.248 
bytes 

No. of bytes in distributed program, including test data, etc: 35.636 bytes 
Keywords: Slavnov-Taylor, BRST, Mathematica 

Nature of physical problem: Symbolic computation in the Slavnov-Taylor for- 
malism for gauge theories in 4-Dimensional space-time based on a semi-simple 
compact Lie group for a general BRS transformations. 

Restrictions on the complexity of the problem: Only matter in the adjoint is 
allowed. 

Typical running time: less than one second 
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2 Introduction 



Symbolic computation in the Slavnov- Taylor (ST) formalism is very useful 
to study properties of a field theory such as quantum stability. The main 
computation problem is to find the most general counter-term compatible 
with the ST identity, and then absorb this counter-term by a field and coupling 
constant redefinition. This can be a not easy task if the BRS transformation 
is not nilpotent (0). It is therefore desirable to construct a calculational tool 
which may perform the computation automatically. In this article we present a 
tool (Slavnov-Taylor ) which partially solves this problem. Slavnov-Taylor 
allows us to perform generic computation of BRST variation and construct 
generic polynomials of a given dimension and ghost number. It can be used to 
check hand computation for generic ST operator which may not be nilpotent. 
Moreover the symbolic manipulation allows us to perform computation in 
general gauges. 

We applied Slavnov-Taylor to verify the correctness of the computation used 
to prove the quantum stability of the Curci-Ferrari gauge in a gluon-ghost 
condensation scenario (0), where the BRS is not nilpotent. 

2.1 The package 

Slavnov-Taylor has been written as a Mathematica package. The interface 
with the user is given by the Mathematica application, which can be run either 
in on-line command or with the front-end feature. 

The algorithm is not optimized, because the running time is not an important 
point due to the fact that people will need to run it few times (or even once) 
for a given model. 

The imposition of the ST identity turns out to be equivalent to the solution of 
a linear problem. Slavnov-Taylor is a tool whose output could be used for 
constructing such linear system. 

It may also be used as help for more theoretical studies such as the stability 
of gauge theory when the ST is not nilpotent. 

2.2 Aims and Contents 

The main aims of this article are to provide a manual for the use of Slav- 
nov-Taylor , to describe the limit of application and the notation used (to 
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allow for user generalization). 



The rest of this paper proceeds as follows: 

• The relevant parameters are presented in sec. 3. 

• The approximations employed are noted in sec. 4. The algorithm of the 
calculation is also outlined. 

• A description of how to use the package is given in sec. 5, including infor- 
mation on the input. More technical informations related to running and 
extending the package are placed in appendices. 

• Sample output from one run is displayed and explained in appendix A. 

• Appendix B gives a more technical example of the package. 

• Finally, in appendix C, a description of the relevant modules and objects 
and their relation to each other is presented. 



3 The creation of a model 

In this section, we introduce the objects and parameters in the Slavnov-Tay- 
lor conventions. Translation to the actual name used in the Slavnov-Taylor 
package are shown in Appendix C. 



3.1 Field content 

The gauge theories in 4-Dimensional space-time considered in the Slavnov-Ta- 
ylor package are based on a semi-simple compact Lie group. In a general 
gauge model the gauge fields are Lie algebra valued and the belong to the 
adjoint representation: 



Matter fields can be included if they are in the adjoint representation (which 
is for example the case of the gaugino for pure Yang-Mills theory in Super- 
symmetry) . 

The action is constructed as the most general gauge invariant and power- 
counting renormalizable action. We define the curvature: 




[7"a,T6] = ifabc'Tc 



C-i 



TrTaU = 5ab 
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and the invariant is given by: 

S- = — f d^T F F^" 

'-'mv J /ii^--' 

The operator . is the trace over the gauge group index. 

The gauge fixing is done in the BRS way: one introduce the Faddeev-Popov 
ghost and antighost fields c(a;)" and c{x)"', and the Lautrup-Nakanishi field 
b{x). The latter normally play the role of a Lagrange multiplier for the gauge 
condition. These three fields are Lie-algebra valued: c{x) — c"(x)t", etc. They 
are used to write the gauge fixing term, needed to extend to quantum level 
the theory. 

Renormalization requires the introduction of an external source coupled to the 
BRS variation of the fields. This is done by adding to the action the terms 

Sext^ [ d'x Y.X[^.{S^) 

where the sum is over the fields which have a non linear BRS transformation. 
We used the Slavnov-Taylor notation as the external field coupled to 

the non linear BRS variation of the $ field. 

3.2 BRST transformation 

The BRS transformation are defined as: 

1 

sc = i-c A c 
2 

sc = b 
sb = 

where the A operator is the external product and is defined on two general 
quantities Lie algebra valued x ^^^^ V 

iX^Vr^fabcXV 

BRST invariance can thus be expressed by a functional identity, the ST iden- 
tity: 

S{S) = 

where the nonlinear ST operator is given, for any functional F by: 
Q^p^ / ,4 V- 5F 
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where the sum over $ is a over the fields with non hncar BRS transformation, 
while the sum over (p is over the fields with linear BRS transformation. 
In case the BRS is nilpotent the gauge fixing action is given by an s variation 
of an integrated local polynomial: 



otherwise it is a ST invariant integrated local polynomial. 
ST formalism allows us to study stability and unitarity of the model, but these 
properties are not investigated by our package. These points are reserved for 
a future version of the Slavnov-Taylor package. 

3.2.1 Fields Properties 

The field properties used by the Slavnov-Taylor package are: 

3.2.1.1 Lorentz properties: The boson gauge fields A^{x) have a Lorentz 
index, while the b"'{x), c{x) and c{x) have no Lorentz index (they are scalar). 

3.2.1.2 Statistic: The boson gauge fields A^(a;) and the Lautrup-Nakanishi 
fields b"'{x) are commuting fields, while the components of c{x) and c{x) are 
anti-commuting. 

3.2.1.3 Canonical dimension: The boson gauge fields ^^(a;), the ghosts 
c"(x) and the antighost c"(x) have canonical dimension one, while the fields 

have canonical dimension two. 

3.2.1.4 Ghost number: The boson gauge fields ^^(a;) and the Lautrup- 
Nakanishi fields are fields with ghost number zero, while the components 
of c{x) have ghost number one and the components of c{x) have ghost number 
minus one. 



4 Calculation 

In fig 4 we show the algorithm used to perform the calculation. Due to the fact 
that user can define the field content he needs and their BRS transformation, 
he can easily modify the model. For example to study pure Supersymmetric 
Yang-Mills theory, he needs to define gaugino fields and the Supersymmetric 



Sgf = s d^x P {fields) 
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User loads the Slavnov-Taylor package 



User defines the Fields 



User defines the BRS transformations 



Package defines the Anti-fields 



Package defines external action 



User obtains ST transformations of a general integrated polynomial 



and the general polynomial of given dimension and ghost number 



Fig. 1. Algorithm used to calculate ST transformations of general quantities in a 
given model. Each step (represented by a box) is detailed in the text. 

BRS transformations. In the same way the generalized BRS transformations 
introduced in (j^ can be treated. Moreover theories with non nilpotent BRS 
transformation can be taken into account with no problem (for example on- 
shell formulation of Supersymmetric theory (0), or gauge theory in presence 
of condensates 



5 Running Slavnov-Taylor 



The program ST.nb is included in the Slavnov-Taylor distribution. This 
program is can be run by Mathematica. The program contains: 

• the definition of the fields a, c, c, b and their BRS transformation; 

• the invariant action and the gauge fixing Curci-Ferrari action as an s vari- 
ation of a polynomial; 

• The definition of the full action as the sum of the invariant action, the 
gauge-fixing action and the external action. 

The program gives as output: 

• the external action Sext create by the program; 

• the list of the field defined by the user; 

• the list of the external field defined by the program; 

• the explicit form of the full action; 
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• the dimension of the full action; 

• the ghost number of the full action; 

• the values of c. 

• the s variation of the full action (which is zero); 

• the most general polynomial of dimension three and ghost number —1; 

• the most general polynomial of dimension two and ghost number 0; 

• the s variation of af^.a^; 

• the s variation of the s variation of a^.a^] 

5.1 Description of the main procedures and global parameters 

Setting the global parameter $ QUIET to a non-zero value gives additional 
information on field creations. 

The properties of the fields content of the model arc encoded in objects and 
must be provided. The same for the BRS transformations. We give here a 
description of the main procedures that user can use. 

• DefineField 

• Parameters: 

Name is the name of the field; 

Indices is a list which determines the Lorentz properties of the 
field (an empty list define a scalar); 

Type determines the statistic of the field, Type=$BOSE define a 
bosonic (commuting) field; Type=$FERMI define a fermionic (anti- 

commuting) field; 

Dimension is the dimension of the field; 
Gh define the ghost number of the field. 

• Objected created: 

Field 

• Set quantities: 

Stat the statistic of the field: +/ — 1 for commuting/anti-commuting 
field; 

Dim the canonical dimension of the field 
Gh the ghost number of the field 

• Global variables modified: 

appends to the variable $FieldList the name of the field 

• SetBRST 

• Parameters: 

Name The field name of which we are defining the BRS transforma- 
tion; 

brs The BRS transformation of the field. 

• Object created: 

$ExtField (if the BRS transformation is not linear in the fields); 
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• Set quantities: 

S is set to be the BRS transformation; 
BRST is set to TRUE; 

• Global variable modified: 

for non linear BRS transformations it adds to Sext the right quantity 

• Module called: 

for non linear BRS transformations it call the module Def ineExtField 
(which appends to the variable $ExtFieldList the name of the 
external field, and works in a similar way than Def ineField mod- 
ule) 
• Create 
■ Parameters: 

Dimension the canonical dimension of the integrated polynomial 
we are looking for 

Gh the ghost number of the integrated polynomial we are looking for 

An object Field can be removed by the module Undef ineField. An 
object ExtField can be deleted by the module Undef ineExtField. 



A Sample Output 

We now present the output for a single example calculation. When you load 
the package you obtain some information on the version in use. We load the 
gauge theory model with usually BRS transformation which is defined in the 
ST . nb file. Once the model is defined, we can perform the computation. We 
run the commands 

In[l]:= «PackageST.m; 
In [2] := «ST.iib; 
In [3] := 

s [a [mu] . a [mu] ] 

s [%] 

The output obtained is 

mu mu 
Out [3]= -2 a . d c 
Out [4]= 

The first output is the consequence of the application of the ST operator to 
quantity a^Ta * T5. The second output is the consequence of the application 
of the ST operator twice to the same quantity. In this case, due to the nil- 
potency of the ST operator, the result is zero. 
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Table B.l 



Detailing arguments in order for module used to define a model. 



Name 


Arguments 


Def ineField 
SetBRST 


Name, Lorentz indices, Statistic, Dimension, Ghost Number 
Name, BRST transformation 



B Sample Program 



We now present the sample program from which it is possible to run Slav- 
nov-Taylor in a simple fashion. The most important features of the modules 
and objects are described in appendix C. The sample program has the form 
displayed in figure B. 

After an initial introductory print-out, the fields a,b,c,c are defined and their 
BRST transformation are set. For these, the same notation as appendix 5 is 
used. 

The user has created the field content of his model. Now is time to define the 
classical action and to fix the gauge (for example the Curci- Ferrari gauge): 
The first results that can be easily obtained are the external action Sg^u the 
list of defined fields and external fields. Moreover the explicit form of the full 
action is obtained, and the dimension and ghost number of a general integrated 
polynomial (in the example Sfot)- To show the use of the ST operator s, an 
application of the Dr operator is shown and the variation of the full action 
(which turn out to be zero) is reported. 

The use of the module Create is shown and the most general polynomials 
of dimension 3 and ghost number -1, and dimension 3 and ghost number zero 
are reported. 

Finally a simple example showing the nilpotency of the s operator is reported. 



C Modules and objects 



In this section we stretch the most important features of the modules and 
objects. In tab. C.l we list the modules with their parameter with a short 
description, in tab. C.2 and tab. C.3 we give a list of the operators with their 
arguments and description. 
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Table C.l 

List of the modules with their parameter and description 



Name 


Arguments 


Output 


Gh 


Pol 


Ghost number of Pol 


UndefineField 


name, Lorentz 


Delete a Field 


DefineField 


name, Lorentz, type, dim, gh 


Create a Field 


UndefineExtField 


name, Lorentz 


Delete an ExternalField 


DefineExtField 


name, Lorentz, type, dim, gh 


Create an ExternalField 


SetBRST 


name, BRS 


Associate a BRS variation to a Field 


Create 


dimax, gh 


Generate the most general integrated poly- 
nomial of the fields with dimension less 
or equal to dimax and ghost number gh 



Table C.2 

List of the operators with their arguments and description 



Name 


Arguments 


Output 


Stat 


Mon 


+1 if Mon is bosonic, -1 if it's fermionic 


Dim 


Pol 


The dimension of Pol 


Dr 


Pol, name 


The functional derivative of Pol with respect to field name 


s 


Pol 


ERST variation of Pol 



Table C.3 

List of the overloaded operators with their arguments and description 



Name 


Arguments 


Output 


NonCommutativeMultiply 

Wedge 
Dot 


Poll, Pol2 

Poll, Pol2 
Poll, Pol2 


Take into account the non commutativity of 
Lie algebra valued polynomial 
External product between Poll and Pol2 
Scalar product between Poll and Pol 2 
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ST.nb 



Automatic computation in tlie BRST formalism 



■ Definitions 

■ Reading Package 

In [1] := « PackageST.m; 

Package ST 

Version 1.0 by Marco picariello 
In [2]:= $Quiet = False; 



■ Definition of tlie fields 

(jq" In [31:= DefineFi6ld[a, {n}, Boss, 1, 0] 

DefineField[b, {}, Bose, 2, 0] 
^ DefineFieldEc {), Fermi, 1, 1] 

• Def ineField[barc, {}, Fermi, 1, -1]; Format[barc] = c; 

Defined field a""- 
tr* Defined field b 

rr t/! Defined field c 
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Defined field bare 



3 

■ Definition of the transformation of tlie fields 

O 

mni:- SetBRST[c, (1/2) cac] 
fO SetBRST[barc, b] 



SetBRSTIb, 0] 
SetBRST[a[n_] , - Alii, c] -a[;u]Ac] 

Defined classical field x[c] 
Defined classical field x[a]"'"- 

■ Definitions of the action 

■ Definition of ttie constants 

In[ll]:= MyNumberQIg] = True; 

HyNumberQ[a] = True; 

■ Definition of the classical action 

in[l31.- F[i_, j_] =d[i, a[J]] -d[j, a[i]] + a[i] Aa[j]; 
Sinv = l/g*2F[/i, v].F[/i, v] ; 

in[l51:- Sgf = s [ bare . d[;i, a[|i] ] + ( a / 2) bare, (b - (1 / 2) bare a c) ] ; 

■ The full action 

In[16]:= Stot = Siil5Jlify[Expand[ (Sinv + Sgf + Sext) ] ] ; 



ST.nb 
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■ Results 

■ The action 

mini:- Sext 

Out [171- (d" c) .X [a]" - c. (a" A x[a]") + y c. (c a x[c] ) 

ln[18]:= $FieldList 
Out[181= {a"'-*, b, c, c} 

In[191 := $ExtFieldList 

Out[191- {X[c] , X[a]"") 

In[201:- StOt 

Out[201- "^"^''>-"^'^^' - ^ (d^a'-) (d-a-^) ^ „ ^ i<i^an.«i-'.n _ 

^""■"y - 2j^^AllL^_U^U^,^. (d. a") + ^ - i ab. (CAC) . 

c. (d" (d'' c) ) + c. (a*' « (d*' c) ) - c . (a" a x [a] " ) - c. (c a (d'' a" ) ) + 

1 , rni (a^.a^)**(a^.a^) (a'^.a^)**(a^.a^) 1 , , , , 
y c. (c aX[c] ) + ^ '-^^ 5 '-^ + — o (c.c) (c.c) 

In[211:- Dim[StOt] 

Out [211= 4 

In[22]:= Gh[Stot] 
Out [221= 

■ Dr operator and how to use it 

In[231:- e.Dr[Stot, b] 

Out [231= a b.c + c. (d^ a'') - y a c. (c a c) 

■ Computations s[Stot] 

ln[241:= Simplify [s [Stot] - ot / 8 c. Jacobi [ (bare a bare) , c, e] ] 

Out [241= 

■ General Counterterm 
in[251:- Create[3, -1] 

Out [251- kt (d™ a"") .0 + k5 a™ . (d"" 0) + k3 a™ .X [a]™ + k6 b. + k2 c.x [c] + k7 c. (c a c) 
in[261:- Create[2, 0] 
Out[261- k9 a™.a"" + klO CO 



■ Example of the nilpotence of s 

ln[271:- s [a [mu] . a [mu] ] 
s[%] 

Out [271- -2 a"" . (d"" c) 
Out [281- 



